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The star v Andromeda is orbited by three known planets (1), the innermost 
of which has an orbital period of 4.617 days and a mass at least 0.69 that of 
Jupiter. This planet is close enough to its host star that the radiation it absorbs 
overwhelms its internal heat losses. Here we present the 24- ^m light curve 
of this system, obtained with the Spitzer Space Telescope. It shows a clear 
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variation in phase with the orbital motion of the innermost planet. This is the 
first demonstration that such planets possess distinct hot substellar (day) and 
cold antistellar (night) faces. 

Last year, two independent groups (0121) reported the first measurements of the infrared light 
emitted by extrasolar planets orbiting close to their parent stars. These 'hot Jupiter' © planets 
have small enough orbits that the energy they absorb from their hosts dominates their own inter- 
nal energy losses. How they absorb and reradiate this energy is fundamental to understanding 
the behavior of their atmospheres. One way to address this question is to monitor the emitted 
flux over the course of an orbit, to see whether the heat is distributed asymmetrically about the 
surface of the planet. 

We have observed the v Andromeda system with the 24- (xm channel of the Multiband Imag- 
ing Photometer for Spitzer (MIPS) © aboard the Spitzer Space Telescope (7). We took 168 
three-second images at each of five epochs spread over 4.46 days (97% of the 4.617-day or- 
bital period of v Andromeda b) beginning on 18 February 2006 at 12:52 UTC (Coordinated 
Universal Time). After rejecting frames with bad pixels near the star and those with Spitzer's 
'first frame effect' (0) (2% - 8% of the data, depending on epoch), we measure the flux of the 
system and that of the surrounding sky using both sub-pixel, interpolated aperture photometry 
and optimal photometry (8 9) on each frame. 

The detection of eclipses UTk from the "hot Jupiter" planetary systems HD 209458b ©, 
TrES- 1 © and HD 1 89733b (11) demonstrate that a small fraction (~ 0.1%) of the total infrared 
light we observe from these systems is actually emitted from the planet rather than the star. 
Thus, if we can measure the flux of a system at a signal-to-noise ratio (S/N) > 1000, temperature 
differences between the day and night faces of the planet will appear as an orbital modulation 
of the total system flux. With a star as bright as v Andromeda, our three-second exposures each 
have S/N ~ 500, so that our SNR expectation is ~ v / 160 x 500 « 6300 at each epoch. 

The MIPS instrument acquires data by placing the stellar image in a sequence of 14 positions 
on the detector. The detector's response varies with position at the ~ 1% level. This variation 
is stable and reproducible, so we calculate correction factors as follows: At each epoch we 
compute the mean measured system flux at each position and take the ratio with the mean in 
the first position. We then average this ratio over all epochs for each position. This results in 
corrections < 2% between positions, with uncertainties ~ 6 x 10~ 4 . Bringing the photometry 
to a common normalization allows us to average over all the frames in each epoch to achieve 
S/N ps 4350 at each epoch. 

As with most infrared instruments, MIPS's sensitivity varies in time. We correct for such 
drifts by dividing the system flux value by the measured background in each frame. The back- 
ground at 24 |xm is thermal emission from the zodiacal dust. This dust pervades the inner solar 
system, absorbing light from the sun and reradiating it at infrared wavelengths. At 24 |xm, its 
emission is strong enough for use as a flux standard, a technique used successfully in measuring 
the eclipse of HD 209458b ©. However, the present work requires one additional correction. 
The zodiacal background is the integrated emission by dust along the line of sight between the 
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telescope and the object. The observed value thus undergoes an annual modulation as that line 
of sight varies with the telescope's orbit about the sun. The best available model (12) predicts 
a linear drift over the brief interval of our observations. However, we cannot use the Spitzer 
model directly, since it is calculated for a line of sight from the Earth to the object in ques- 
tion. The difference in position between the Earth-trailing telescope and the Earth itself is large 
enough that the slope of the variation may be slightly different. Thus, we fit for the linear drift 
directly, simultaneously with any model lightcurve fits. 

The phase curve for the v Andromeda system shows a variation (Fig [l]) in absolute pho- 
tometry, even before any corrections for instrumental or zodiacal drifts are made. After the 
calibration with respect to the zodiacal background is applied, this variation is revealed to be in 
phase with the known orbit of the innermost planet of the system, our principal result. 

A simple model can be fit to the phase curve (FigO, assuming local, instantaneous thermal 
reradiation of the absorbed stellar flux. In the simplest model, the phase of the variation is not a 
free parameter, but is rather set by the measured radial velocity curve (2). Phase offsets are pos- 
sible for models in which the energy is absorbed deep within the atmosphere and redistributed 
about the surface U3\\14\i . There is weak (2.5a) evidence for a small phase offset in this data 
(Fig El) but the large offsets predicted from some models are excluded at high significance. Fit- 
ting the peak-to-trough amplitude to the observations yields a best-fit value for the planet-star 
flux ratio 2.9 ± 0.7 x 10~ 3 times the star's brightness. This is very similar to the result at this 
wavelength for HD 209458b (4). However, the latter is a measure of the absolute flux from 
the planet divided by that from the host star, while the present result is a measure of the flux 
difference between the projected day and night sides, divided by the flux of the (different) host 
star. 

Another difference between the cases of v Andromeda b and HD 209458b is that we do 
not have a strong constraint on the orbital inclination in this system, so we must include the 
unknown inclination in the model fit (Fig O • At higher inclinations, parts of both the night side 
and the day side are always visible, so the true contrast between the day and night sides must 
be larger than the amplitude of the observed variation. This contrast is ultimately driven by the 
light absorbed from the star, which therefore provides an upper limit. We know the distance of 
the planet from the star and the stellar properties, so we can estimate the contrast that would 
result if all of the observed flux were reradiated from the day side and nothing from the night 
side. If we assume the planet's radius is <1.4 Jupiter radii (as observed for other planets of this 
class), then we can constrain the expected amplitude to be < 3.4 x 10 -3 (2a) for a simple black- 
body, no-redistribution model with zero albedo. Thus, a consistent picture of the atmospheric 
energetics emerges as long as the orbital inclination is >30°. 

A natural question to ask is whether there are any plausible alternative models for the ob- 
served variation. The estimated rotation period of the star is too long to explain our phase curve 
as the result of a normal starspot (which is darker than other parts of the stellar surface). One 
could posit a feature on the stellar surface similar to a starspot but induced by a magnetic in- 
teraction between the star and the planet, and therefore moving synchronously with the planet. 
However, (T73T) place an upper limit of 1.6 x 10 -4 on the amplitude of optical variation with the 
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planetary orbital period, so infrared variability from the star should be even weaker than this. 
Some evidence for such magneto spheric interactions is found in observations of chromospheric 
calcium H and K lines (173b and has even been seen in the v Andromeda system. However, the 
energy input needed to explain the Ca lines is ~ 10 27 ergss _1 , much less than the minimum 
planetary luminosity we infer here (~ 4 x 10 29 ergs s _1 ). Indeed, one can make a quite general 
argument that our observations cannot be powered by the same mechanism, since any heating of 
the star due to magnetic interaction with the planet ultimately extracts energy from the planetary 
orbit. Thus, one may calculate an orbital decay time 
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where M* and M p are the stellar and planetary masses, a is the semi-major axis and E is the 
observed heating rate. Heating at the level necessary to explain our observations would result in 
the decay of the planetary orbit on timescales < 10 7 years, while the estimated age of the system 
is 3 Gyr. As such, the chromospheric heating of the star is unlikely to be related to the effect 
seen at 24 \im. 

This observation reveals the presence of a temperature asymmetry on the surface of an 
extrasolar planet. The first measurements of eclipses (14151) yielded measurements of the absolute 
flux levels emerging from the day sides of two extrasolar planets. When compared with models 
of radiative transfer in such atmospheres lZ7lli#ll7Pl 12(71) . those observations are consistent with 
a situation intermediate between no redistribution and full redistribution. A similar comparison 
is possible in this case (Fig|4|). Our observed day-night flux difference is comparable to the flux 
emerging at full phase in the models of dT7b . which indicates that there is little evidence for 
redistribution of energy to the night side. 

In conclusion, the observation of the phase curve of v Andromeda b indicates that signifi- 
cant temperature differences exist between the day and night faces of the planet, consistent with 
a model in which very little horizontal energy transport occurs in the planetary atmosphere. 
Furthermore, it indicates that the opportunities for direct extrasolar planetary observations are 
better than previously thought, since useful data can be obtained even in cases where the plan- 
etary orbit is not so fortuitously aligned that the system exhibits transits or eclipses. 
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Figure 1 : The light curve of the v Andromeda system. A) The phase variation in the v An- 
dromeda system flux before any corrections are applied for instrument or zodiacal drifts. Vari- 
ations in the system flux are significant even at this point. B) By comparing to the zodiacal 
background, and fitting for the linear drift in the background due to the telescope's motion, we 
obtain the phase curve shown. In each case, phase is shown modulo unity, with zero phase 
occurring when the planet is closest to earth. The amplitude units are expressed in terms of the 
system flux at the first epoch. 
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Figure 2: Comparison of the phase curve and the No-Redistribution Model. The solid points 
show our final phase curve, after applying calibrations, in time order from left to right. The open 
points are repetitions of these, displaced horizontally by one orbit, to better illustrate the phase 
coverage over two cycles. The solid line is an analytic model for the planetary emission in which 
energy absorbed from the star is reradiated locally on the day side with no heat transfer across 
the surface of the planet, the so-called No-Redistribution model (and in excellent agreement 
with the more detailed version in (I7%1) ). The assumed inclination in this case is 80° from pole- 
on, and the relative planet/star amplitude is 2.9 x 1CT 3 . If we allow for a phase shift relative to 
the radial velocity curve, we obtain a slightly better fit, as shown by the dotted curve. The best 
fit is obtained with a phase lag of 11°, but zero lag is excluded only at the 2.5a level. 
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Figure 3: The influence of inclination on the inferred Day-Night contrast. The solid contours 
bound the 1, 2 and 3a confidence regions for the day-night flux difference (in units of the stellar 
flux), determined as a function of assumed orbital inclination (measured relative to a face-on 
orbit). The large shaded regions indicate those values excluded at 3 a. The lower shaded region 
is excluded because the planet does not transit in front of the star. The vertical dashed line indi- 
cates the expected upper limit to the contrast, obtained when the night side is completely dark 
and all of the stellar flux is reradiated from the day side, in accordance with the no-redistribution 
model, and assuming zero albedo. At the right we show the true mass of the planet given the 
assumed inclination (based on the minimum mass derived from the radial velocity curve), in 
units of Jupiter masses. 
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Figure 4: Comparison of the measured amplitude and a planetary spectral model. A) The 
solid curve shown is a model (f77b for a planet of radius 1.4 irradiated with parameters 
appropriate to the v Andromeda system observed at full phase. This results in a temperature 
~ 1977 K (TZ7I). The model is in agreement with the observations (filled circle) at the 2a level. 
B) The normalised spectral response curve of the MIPS 24 fim instrument extends from 20fim 
to 30/um. 
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